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SECTION  i 


INTRODUCTION 

In  many  cases,  it  is  helpful  to  calculate  the  failure  probability 
for  an  aircraft  in  order  to  obtain  a  more  quantified,  indication  of  the 
risk  involved  with  a  particular  course  of  act Those  types  of 
calculations  have  proven  useful  in  the  past  for  helping  managers  to 
select  from  some  difficult  options. 

The  process  defined  in  this  report  presumes  that  the  analyst  has 
the  following: 

(1)  The  probability  density  function  for  crack  length  for  an 
area(s)  of  the  aircraft  that  is  to  be  analyzed. 

(2)  The  probability  density  function  for  stress  for  the  area(s) 
of  the  aircraft  that  is  to  be  analyzed. 

(3)  The  stress  in  the  analysis  area(s)  at  which  a  given  crack 
size  would  be  critical. 

(4)  The  inspection  raliability  for  inspections  that  are  to  be 

made . 

(5)  The  structural  rework  to  be  performed  (which  consists  of 
fastener  hole  reaming). 

(6)  The  crack  growth  in  th>  analysis  area(s)  as  a  function  of 

time . 

(7)  The  mission  time  for  the  desired  sequence  of  flights. 

This  information  is  input  to  the  computer  routine  as  detailed  in 

Section  III.  The  output  of  the  computation  is  the  single  flight, 
probability  of  failure,  the  aircraft  failure  probability,  and  the  expected 
number  of  failures  in  the  force.  The  detailed  equations  for  the 
computations  are  given  in  Section  II. 
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Generally,  the  single  flight  probability  of  failure  and  the 
expected  number  of  losses  are  the  more  useful  outputs  of  this  analysis. 
For  an  in-depth  discussion  of  acceptable  lovels  of  safety,  the  reader 
is  referred  to  Bo  Lundberg's  Wright  Brothers  Lecture  on  "Fatigue  Life 
of  Airplane  Structures"  published  in  the  Journal  of  the  Aeronautical 
Sciences,  June,  1955.  Lundberg  provides  some  arguments  for  safety 
levels  that  could  be  used  for  commercial  aircraft.  The  maximum  allowable 
single  flight  probability  of  failure  for  a  military  aircraft  is 
dependent  on  such  factors  as  the  mission  of  the  aircraft,  safety  of 
the  air  crew,  and  feasibility  of  modifications  and/or  inspections  to 
reduce  the  risk. 


SECTION  II 


DEVELOPMENT  OF  HISK  ANALYSIS  METHODOLOGY 

1  BASIC  FUNCTION  DEFINITIONS 

Suppose  that  T  is  a  positive  number  and  that  Np  is  the  number 
of  control  points  on  the  airframe,  N^  is  the  number  of  different 
missions  flown,  Nj.  is  the  number  of  inspections  in  the  time  interval 
[0 , T] ,  NR  is  the  number  of  structural  reworks  to  be  performed  in  the 
time  interval  [0,T],  and  each  of  b,  c,  d,  and  e  is  a  positive  integer 
such  that  b  is  in  [l,Np],  o  is  in  [l,N^],  d  is  in  [0,Nj],  arid  c  is  in 
[0,NR].  Further,  suppose  that  pb$e  is  a  simple  surface  such  that  the 
point  (x,t,  pb$e(x,t))  is  a  point  of  pb$e  only  if  x  is  a  positive 
number,  t  is  a  positive  number  £,  T,  and  pb£e  (x,t)  is  tho  crack 
length  probability  density  between  the  dth  and  d+1  inspection  and 
between  the  eth  and  e+1  rework  (note  that  when  d  =  0,  the  aircraft 
has  not  received  an  inspection  and  when  e  =  0  the  aircraft  las  not 
been  reworked)  for  tho  crack  length  x  at  the  bth  control  point  of 
the  aircraft  for  the  flight  time  t.  pb$e  is  called  the  crack  length 
probability  density  function  corresponding  to  the  bth  control  point, 
dth  inspection,  and  eth  rework. 

*  bede 

Now,  suppose  that  a  is  a  simple  graph  such  that  the  point 

,,  bede,,..  ,  .  ,  ,,  bode  ,  ,  „  ,  „  bede,  .  ,  , 

(t,  a  (t))  is  a  point  ox  a  only  if  each  of  a  (o)  and  t  is 

bCQG 

a  positive  number  and  a  (t)  is  the  crack  length  at  the  time  t 

bed  c 

(based  on  an  initial  crack  of  size  a  ( o ) )  at  the  bth  control  point 
of  an  aircraft  flying  the  eth  mission  at  the  time  t  after  the  dth 
inspection  and  eth  rework.  The  simple  graph  aboc*e  j.8  called  tho 
crack  growth  function  for  the  bth  control  point  and  the  eth  mission 
and  for  tho  dth  inspection  and  eth  rework. 

Now,  suppose  that  PTb  is  a  simple  graph  such  that  the  point 

U  Jb  Jj 

(x,  Pj  (x))  belongs  to  P^  only  if  x  is  a  positive  number  and  Pj  (x) 
is  the  probability  that  a  crack  of  length  x  will  not  be  detected 
during  an  inspection  of  tho  bth  control  point.  The  simple  graph  Pjb 
is  called  the  inspection  function  For  the  bth  control  point. 
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Suppose  that  p^  is  a  simple  graph  such  that  the  point  (x, 

pgu(x))  belongs  to  pgC(x)  only  if  x  is  a  number  and  pgC(x)  is  the 

probability  density  of  the  bth  control  point  for  the  maximum 

be 

stress  x  in  a  single  flight  of  the  cth  mission,  pg  is  the  stress 
probability  density  function  for  the  bth  control  point  of  an  aircraft 
flying  the  cth  mission. 

b 

Next,  suppose  that  s“  is  a  simple  graph  such  that  the  point 
b  b  c  b 

(x  S°(x))  belongs  to  Sc  only  if  Sc  is  the  stress  at  which  a  crack  of 
length  x  at  the  bth  control  point  becomes  critical.  is  the 
critical  oraok  length  funotion  for  the  bth  control  point. 

Now,  suppose  that  Np  is  a  simple  graph  such  that  the  point 
(t,  Np( t ) )  belongs  to  Np  only  if  t  is  in  [0,T]  and  Np(t)  is  the 
number  of  flights  flown  at  the  time  t.  Np  is  called  the  flight  number 
function. 

Further,  suppose  that  NMp  is  a  simple  graph  such  that  the  point 
(F,  ( f ) )  belongs  to  NMp  only  if  f  is  a  positive  integer  in  the 

interval  [1,  Np(T)]  and  Nwp(f)  is  'the  mission  number  flown  on  the 
flight  numbered  f.  N^p  is  called  the  mission  number  function. 

Suppose  that  tp  is  a  simple  graph  such  that  the  point  (f,tp(f) 
belongs  to  tp  only  if  g  is  an  integer  in  [l,Np(T)]  and  tp(f)  is  the 
time  just  preceding  the  fth  flight.  tp  is  called  the  flight  time 
funotion. 

Also,  suppose  that  Npw  is  a  simple  graph  such  that  the  point 
(b,Npw(b))  belongs  to  Npw  only  if  b  is  in  [l,Np]  and  Npw(b)  is  the 
number  of  points  with  the  same  p^e,  a*30,  p^c,  and  functions 
in  the  computation  for  the  bth  control  point  Np,^  is  the  repeat  control 
point  function. 
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2  INFLUENCE.  OF  INSPECTIONS  ON  THE  CHACK  PKOBAB1I.1TY  DENSITY  FUNCTION 


3  INFLUENCE  OF  STRUCTURAL  REWORK  ON  THE  CRACK  PROBABILITY  DENSITY 
FUNCTION 

Suppose  that  each  of  aRRo  and  aRR  is  a  positive  number  such  that 

0  ^  aRR  —  aRRo‘  T,ie  craci<  length  probability  density  function 

corresponding  to  the  bth  control,  dth  inspection,  and  nth  structural 

rework  is  modified  by  structural  rework  as  follows:  For  each  point 

(x,y )  in  the  x-y  projection  of  p^^such  that  if  x  >  then  x  is  replaced 

b  L  —  HRD 

by  x-aRR.  Therefore,  if  after  this  replacement 


Y  -  / 


PRo  bdn ,  b 
P  r  <aot 


•t)dlx 


then , 


(Ix,t)  dlx 


bdn+i,  ,  . 
p  r  (x,t) 


bdn.  . . 
P  C  (aRRo’tJ 


for  x  <.  a!i 


bdn+1,  , .  bdn,  ,  . 
p  c  (x,t)  =  P  c  (x,t) 


for  x  ^  aRRo  , 


It  is  noted  that  the  new  distribution  is  uniform  between  0  and 
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4  CALCINATION  OK  IMF,  PKOBAh'lI.l  i’Y  OK  FAll.URL 


At  a  given  tim*  the  single  flight  probability  of  failure  can 

be  computed  as  follows:  If  it  is  assumed  that  at  a  given  time  the 

aircraft  crack  distribution  function  is  statistically  independent 

of  the  stress  density  function  for  a  particular  mtrol  point  then  if 

t  is  in  [0,T] ,  b  is  in  [l,Np],  c  is  in  [1,N^],  is  in  [ 1 .  ] »  e  is 

be 

in  [1,N0],  s  is  in  the  x-projection  of  p  ,  x  is  in  the  x-projection 

bde*'  ^ 

of  p  r  then  there  is  a  number 


bdee.  , .  bde,  ,  ,  be 
Pj  (x.s.t)  -  p  c  (x,t)  •  ps  (g) 

which  is  the  joint  probability  density  of  crack  length  and  stress 
for  the  point  (x,s,t).  Since  this  produce  exists  for  each  number 

in  [0,1],  in  the  x-projection  of  pb^e  and  in  the  x-projection  of  pbc 

bede  L  ° 

then  a  simple  surface  Pj  called  the  j.  int  probability  density 

function  of  crack  length  and  stress  is  defined. 

It  follows  then  that  if  Rb  is  a  point  set  such  that  ordered  pair 

(x,s)  is  a  member  of  Rb  only  if  the  crack  length  is  critical  for  the 

stress  s  (defined  by  the  simple  graph  S  )  then  the  probability  of 

failure  a.  the  bth  control  point  at  the  time  t  during  one  flight  of 

the  cth  mission  is 


bede .  . 

Pv  = 


u 


p^cde[ix,ls,t]dlxdls 


The  probability  that  there  will  be  no  failure  at  the  time  t  for 

bede 

a  single  flight  of  the  cth  mission  is  (l-pp  (t)).  Therefore,  for 
all  control  points  identical  to  the  bth,  the  single  flight  probability 
of  failure  is 

bode 


FW 


...  ,  NPW(b)  bede.,.. 

(t)  =1-77-  U-PF  (t)) 

1 


and  for  the  aircraft  the  single  flight  probability  of  failure  is 


ode 

VWA 


■l-^p  n-Jfttu 


b=l 
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SECTION  III 


DESCRIPTION  OF  COMPUTER  PROGRAM 


1  NOTATION 

The  notation  for  input  data  to  the  routine  is  given  later  in  this 
section  and  is  not  repeated  here. 

TF  Table  -  a  set  of  numbers  used  for  single  linear  interpolation 
to  determine  the  flight  time  corresponding  to  the  number  of  missions 
flown  by  the  aircraft.  The  maximum  number  of  points  allowed  in  the 
TF  Table  is  100. 

PCB  Table  -  a  set  numbers  that  contains  points  of  the  crack 
length  probability  density  function  for  interpolation  in  the  PCBM 
Tabie.  The  maximum  number  of  points  allowed  in  the  PCB  Table  is  51. 

PI  Table  -  a  set  of  numbers  used  for  single  linear  interpolation 
to  determine  the  probability  of  not  detecting  a  crack  of  a  given 
size  during  an  inspection.  The  maximum  number  of  points  allowed  in 
the  PI  Table  is  100. 

SCB  Table  -  a  set  of  numbers  used  for  single  linear  interpolation 
to  determine  the  stress  at  which  a  given  crack  siz..'  will  be  critical. 

The  maximum  number  of  points  allowed  in  the  SCB  Table  is  100. 

PCBM  Table  -  a  set  of  numbers  computed  from  the  PCB  Table  used 
for  single  linear  interpolation  to  determine  the  crack  length 
probability  density  function  for  a  given  crack  size.  The  PCBM  Table 
is  modified  by  the  routine  to  account  for  crack  length  change. 

PDSBC  Table  -  a  set  of  numbers  used  for  double  linear  interpolation 
to  determine  the  stress  probability  density  function  for  a  given  stress. 
Total  number  of  PDSBC  Table  entries  allowed  is  1000.  Maximum  number 
of  stress  entries  per  mission  is  50. 
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AB  Table  -  a  set  of  numbers  used  for  double  linear  interpolation 
to  determine  the  crack  length  for  a  given  mission  and  time.  Used 
to  calculate  the  incremental  crack  growth  for  flight  of  a  given 
mission(s).  Total  number  of  AB  Table  entries  allowed  is  1000. 

Maximum  number  of  craok  length  entries  is  100. 

CRACKL(I)  -  The  Ith  crack  length  entry  in  the  PCBM  Table. 

PCB'i'(l)  -  The  Ith  crack  probability  density  entry  in  the  PCBM 
Table. 

NCRL  -  The  number  of  craok  length  entries  in  the  PCBM  Table. 

DAREAP ( I )  -  Trapezoidal  rule  approximation  to  the  area 

/CRACKL  ( 1+1 ) 

pBC®(Ix,0)dIx 

CRACKL(I) 

PSIL(I)  -  The  Ith  stress  entry  in  the  PDSBC  Table. 

PDPSl(I)  -  The  Ith  stress  probability  density  funotion  entry 
in  the  PDSBC  Table. 

NPSIL  -  The  number  of  stress  entries  in  the  PDSBC  Table. 

TIMCLC(I)  -  The  time  derived  from  the  TF  Table  corresponding 
to  the  number  of  flights  NFLTS(I). 

PFBCT(KT)  -  Probability  of  failure  during  one  flight  at  the  KTth 
risk  calculation  (including  inspections  and  reworks)  for  a  single 
control  point. 

PFBC  =  PFlicr(KT)  for  a  given  KT 
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¥  f . 

TIMI(JINSP)  -  The  time  derived  from  the  TF  Table  corresponding 
■  |  to  the  number  of  f  lights  NIN(JINSP) . 


TIMR(JREWK)  -  The  time  derived  from  the  TF  Table  corresponding 
to  the  number  of  flights  NRE(JREWK). 

TIME(KT)  -  The  time  corresponding  to  the  KTth  riBk  calculation 
(including  inspections  and  reworks). 

DELTM  -  The  incremental  time  between  risk  calculation  points. 

TB  -  The  time  .derived  from  the  AB  Table  corresponding  to  a 
specified  crack  length  and  mission.  The  crack  length  is  an  entry  in 
the  PCBM  Table. 

DLAREA  -  Integration  step  size  for  calculation  of  single  flight 
probability  of  failure. 

PFT(I)  -  The  single  flight  probability  of  failure  for  the 
airoraft  at  the  Ith  risk  calculation  point. 

PF ( I )  -  The  probability  of  failure  for  the  aircraft  for  the  time 
defined  by  the  Ith  risk  calculation  point. 

PFF(I)  -  The  expected  number  of  failures  in  a  force  of  size 
NA  corresponding  to  the  Ith  risk  calculation  point. 


2  COMPUTATION  OF  AIRCRAFT  FAILURE  PROBABILITY 


The  subroutine  FINCLC  computes  the  aircraft  failure  probability 
as  follows:  Suppose  that 

I  is  an  integer 

PFT(I)  is  the  single  flight  probability  of  failure  at  the 
Ith  calculation  step 

PF { I )  is  the  airplane  probability  of  failure  at  the  Ith  cal¬ 
culation  step 

RNFj  is  the  number  of  flights  between  the  I-ist  and  Ith 
oaloulation  steps 

Suppose  further  that  the  single  flight  probability  of  failure 
for  the  RJth  flight  after  the  I-lst  oaloulation  step  can  be 
expressed  by 

PFT];(J)  =  PFT(I-1)  + 

/  PFT(I)  -  PFT(I-lA  rj 
\  RNFj  ' 

Consequently,  the  airplane  failure  probability  at  the  ith 
oaloulation  step  is 

PF  ( I )  =  1.0  -  jp  (  1  " 

RJ  =  1 

(«<>-»  — )  „)) 
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3  COMPUTKK  FLOW  DIAGRAM  AND  PROGRAM 


The  computer  routine  waa  coded  in  FORTRAN  Extended  Language 
with  the  main  program  and  subroutines  arranged  as  follows! 


RISKY  Program  (See  Appendix  A) 


. . 
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MAIN  -  Main  Program  -  Sets  NZERO  to  zero  and  transfers  program 
control  to  GUIDE. 

GUIDE  -  Subroutine  -  Initially  zeros  input  and  output  numbers. 
GUIDE  is  the  primary  controlling  subroutine  that  for  each  of  NP 
oontrol  points  on  the  structure  transfers  control  to  INPUT,  CALC, 
and  PR1NTR. 

INPUT  -  Subroutine  -  Reads  in  all  input  data.  A  procedure  for 
reading  in  the  required  input  data  is  given  later  in  the. section, 

INAID  -  Subroutine  -  Called  by  INPUT  and  has  the  purpose  of 
writing  out  certain  input  data.  INAID  also  sets  NZERO  =  1  for  oontrol 
of  data  handling  in  GUIDE. 

PCBM  -  Subroutine  -  Called  by  INAID  to  modify  the  PCB  table  and 
areate  the  PCBM  table.  If  there  are  n  points  in  the  PCB  table  and 
if  x^,  x2  '  xn  are  the  abscissas  and  y^  y^,  yn  are  the 

ordinates  of  the  points  of  the  PCB  table  then  Xj,  x2,  l.OOOOOlx,,, 
x3,  1.000001x3,  ' 1.000001xn_1,  xn  are  the  abscissas  and 

yl*  y2’  y2’  y3’  y3’  '"yn-r  yn  are  the  ordinates  of  the  Pointa  of 
the  PCBM  table. 

PAGEHD  -  Subroutine  -  Writes  out  page  heading  including  run 
identification,  date  and  page  number. 

CALC  -  Subroutine  -  Performs  the  following  functions: 

(1)  Computes  flight  times  for  failure  probability  calculation. 

(2)  Corrects  the  crack  length  probability  density  function 
for  error  in  area. 

(3)  Identifies  the  mission  being  flown. 
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(4)  Checks  for  end  of  computer  run. 


(5)  Checks  for  inspection  requirement  before  next  flight. 

(6)  Checks  for  rework  requirement  before  the  next  flight. 


(7)  Establishes  new  crack  probability  density  function 
for  next  failure  probability  calculation. 


PCALC  -  Subroutine  -  Integrates  the  joint  probability  density 
funotion  p^1 2 3^0  to  obtain  p'5^6. 


INSPEC  -  Subroutine  -  This  subroutine  is  called  when  an  inspection 
is  to  be  performed.  INSPEC  computes  the  new  craok  length  probability 
density  funotion  after  an  inspection.  Also,  it  corrects  the  orack 
length  probability  density  funotion  for  error  in  area.  It  is  assumed 
that  the  oraoks  found  by  inspection  are  repaired. 

REWORK  -  Subroutine  -  This  subroutine  is  called  when  a  structural 
rework  is  to  be  performed.  REWORK  computes  the  new  orack  length 
density  funotion  after  the  rework  has  been  performed.  Also,  it  oorrects 
the  oraok  length  probability  density  function  for  error  in  area.  The 
rework  performed  is  independent  of  orack  size.  The  crack  population 
is  diminished  by  the  number  aRR  <  aRR0 

FINCLC  -  Subroutine  -  Called  by  CALC  to  perform  the  following 
functions: 

(1)  Computes  the  single  flight  probability  of  failure  for 
the  aircraft. 

(2)  Computes  the  probability  of  failure  for  the  aircraft. 

(3)  Computes  the  expected  number  of  failures  in  the  population. 
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TABS  -  Subroutine  -  Used  for  straight  line  interpolation  in 

a  single  table  look-up. 

TABDS  -  Subroutine  -  Used  for  straight  line  interpolation  in  a 

double  table  look-up.  Suppose  eaoh  of  Sx,  S  ,  and  Sz  is  a  number  set 

such  that  if  (x,y,z)  is  a  point  of  the  double  table  then  x  belongs  to 

Sv,  y  belongs  to  S  and  z  belongs  to  S  .  Further,  suppose  that  S 
a  y  z  xy 

is  a  point  set  in  a  plane  suoh  that  if  x  belongs  to  Sx  and  y  belongs 

to  S  then  (x,y)  belongs  to  S  .  The  interpolation  is  on  S  for  a 
y  _  m  xy  z 

point  (x,y)  that  is  in  the  plane  that  oontains  S  . 

xy 

TABDSI  -  Subroutine  -  Used  for  inverse  interpolation  of  in  a 
double  table  look-up  (see  TABDS).  Suppose  that  S  is  a  point  set 
in  a  plane  suoh  that  if  y  belongs  to  and  z  belongs  to  then 
(y,z)  belongs  to  S  .  If  y  is  in  S  and  z  is  a  number  then  the  inter- 

yz  y 

polation  is  on  Sx  for  a  point  (y,z)  that  is  in  the  plane  that  contains 

syz' 

AREAC  -  Subroutine  -  Computes  area  under  a  simple  graph  based 
on  the  trapezoidal  rule.  Maximum  number  of  points  from  simple 
graph  is  100. 

AREAS  -  Subroutine  -  Computes  area  under  a  simple  graph  based 
on  the  trapazoidal  rule.  Maximum  number  of  points  from  simple  graph 
is  50. 

PRINTR  -  Subroutine  -  Called  from  GUIDE  to  print  the  computed 
results  at  the  conclusion  of  the  run. 


ilftlAVVU.V,  A . WJLM’,1 


M'W.i.i,. 
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4  INPUT  DATA 


All  of  the  input  data  described  below  is  read  into  the  program 
by  means  of  the  subroutine  INPUT.  This  program  is  a  general  purpose 
routine  to  read  the  P  (parameters),  N  (integers) ,  and  the  tables. 
There  are  several  options  by  which  this  may  be  done  by  this  routine. 
A  suggested  arrangement  is  given  as  follows: 


1415  Format 


IDENT 

Li. 

LI 

NPF3 

L°J 

NTH 

NTW1 

MONTH 

DAY 

YEAR 

NTI2 

NTW2 

NA 

NCONTP 

BI5 

Format 

NP  NTMS  | 

NINSP 

|  NREWK | 

;  NMISS  ] 

NLEN  ] 

j  NJST  ] 

|_nprnt] 

72H 

Run ^Description 


72H 


Run  Description 


315  Format 


1 

i  6 

n 

6E10.3  Format 

)  DELCRK  1 

DELST  j 

TIMAX 

AI  ] 

LARR  1 

ARRO 
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If  NTW1  =  0  go  to  (b);  if  NTW1  >  0  go  to  (a) 


Integer  Set  1 


(a) 


315  Format 
| 51  |  55  |  1 


515  Format 


NWT1(1)  through  NWTI(S 


(b)  If  NTW2  »  0  go  to  (d)i  if  NTW2  >  0  go  to  (o) 


Integer  Set  2 


(c) 


31 S  Format 


57  58  1 


215  Format 
NWT2(1)  through  NWT2T2)  1 


(d)  If  NCONTP  <  NP  go  to  (f) 

If  NCONTP  =  NP  go  to  (e) 

(e)  Integer  Set  3 


315  Format 

75  74  +  NP  1 

[  NP  <  15 

1215  Format 

NPW(l)  through  NPW(NP) 
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Integer  Set  4 


Integer  Set  5 


If  NINSP  =  0  go  to  (h) 
If  NINSP  >0  go  to  (g) 


Integer  Set  6 


(h)  If  NREWK  =  0  go  to  (j) 
If  NREWK  >  0  go  to  (i) 


PCB  Table 


15,115  Format 

NSTAB2  <  51  -  NINSP 


6E10.3  Format 

TABLI  (1,2)  through  TABLI  (NSTAB2.2) 
(crack  length  entries) 


6E10.3  Format 

TABLI  (NSTAB2+1 , 2 )  through  TABLI  (2NSTAB2.2) 
(crack  probability  density  entries 


NSTAB2 


PI  Table 


15,115  Format 
Y'  ^  j  N ST  AB 3  ~ 


N STAB 3 <  100 


6E10.3  Format 


TAELI (1,3)  throught  TABLI  (NSTAB3.3) 
(crack  length  entries)  _ 


6E10.3  Format 


TABLI  (NSTAB3+1 ,3 )  through  TABLI  (2NSTAB3.3) 
(inspection  probability  entries) 


3) 


SCB  Table 

15,115  Format 

I  4  )  NSTAB4  ~|  NSTAB4  100 

6E10.3  Format 

TABLI (1,4)  through  TABLI  (NSTAB4 , 4 ) 

(crack  length  entries) _ 

6E10.3  Format _ _ _______ 

TABLI  (NSTAB4+1 ,4)  through  TABLI  (2NSTA34.4) 
(critical  stress  entries) 


(1)  If  NTI2  =  0  go  to  (n) 
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The  first  card  contains  14  fixed  point  (integer)  numbers 
arranged  in  15  fields,  These  14  entries  in  order  on  this  card 
are: 

(1)  1DENT  -  Run  number 

(2)  1 

(3)  0 

(4)  NPF3  -  Number  of  integer  sets  to  be  read  in  -  maximum  number 

(5)  0 

(6)  NTH  -  Number  of  single  tables  to  be  read  in  -  maximum  number 

(7)  NTW1  =  1  if  single  tables  are  to  be  printed 

=  Q  otherwise 

(8)  Month  -  Month  in  date  for  page  heading 

(9)  Day  -  Day  in  date  for  page  heading 

(10)  Year  -  Year  in  date  for  page  heading 

(11)  NTI2  -  Number  of  double  tables  to  be  read  in  -  maximum 
number  =  2 

(12)  NTW2  =  1  if  double  tables  are  to  be  printed 

=0  otherwise 

(13)  NA  -  Number  of  aircraft  in  the  population 

(14)  NCONTP  -  Control  point  number  for  this  run 
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The  second  card  contains  eight  fixed  point  (integer) 
numbers  arranged  in  15  fields,  In  order  these  entries  are: 

(1)  NP  -  Number  of  control  points  on  aircraft  for  analysis  - 
maximum  number  =15 

(2)  NTMS  -  Number  of  times  the  risk  is  calculated  in  the  run 
excluding  the  inspection  and  rework  period 

(3)  NINSP  -  Number  of  inspections  for  this  run 

(4)  NREWK  -  Number  of  reworks  for  this  run 

(5)  NMISS  -  Number  of  different  missions  flown  during  this  run 

(6)  NLEN  -  Number  of  crack  lengths  used  in  failure  probability 
calculation 

(7)  NJST  -  Number  of  stresses  used  in  failure  probability  calculation 

(8)  NPRNT  =  1  if  crack  length  history  is  to  be  printed 

=  0  otherwise 

The  third  and  fourth  cards  each  contain  a  72H  field  for  purpose 
of  run  descriptions,  etc. 

The  fifth  card  contains  1,  6,  and  1  in  15  fields 

The  sixth  card  contains  six  floating  point  number  arranged 
in  E10.3  fields.  These  six  numbers  are  placed  in  the  following  order: 

(1)  DELCRK  -  The  crack  length  interval  for  integration  of  the 

bode 

joint  probability  density  function  pj 
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(2)  DELST  -  The  stress  interval  for  integration  of  the  joint 
probability  density  function  pjjoc*e 

(3)  TIMAX  -  The  flight  time  from  zero  at  the  end  of  the  run 

(4)  AI  - 


(5)  ARR  -  a' 


RR 


(6)  ARRO 


RR° 


The  next  oards  read  in  the  Integer  Sets  as  required. 

Integer  Set  1  - 

First  card  contains  SI,  55,  and  1  in  15  fields 
Second  card  oontains  NWTl(l)  through  NWT1(5)  in  15  fields 
NWTl(i)  a  i  if  the  ith  single  table  is  to  be  printed 
NWTl(i)  s=  0  otherwise 

Integer  Set  2  - 

First  card  contains  57,  68,  and  1  in  15  fields 

Second  card  oontains  NWT2(1)  and  NWT2(2)  in  15  fields  where 
NWT2( i )  =  1  if  the  ith  double  table  is  to  be  printed 
NWTl(i)  =  0  otherwise 

Integer  Set  3  - 

First  oard  oontains  75,  74  +  NP,  and  1  in  15  fields 

The  next  card(s)  contain  NPW(l)  through  NPW(NP)  in  15  fields 

12  entries  per  card  where  NPW(b)  =  Npw(b) 
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Integer  Set  4  - 

First  card  contains  101,  100  +  NTMS,  and  1  in  15  fields 
The  next  card(s)  contain  NFLTS(l)  through  NFLTS(NTMS)  in  15  fields 
12  entries  per  card  where  NFLTS(i)  is  the  number  of  flights  at 
the  ith  risk  calculation  (excluding  inspections  and  reworks)  in 
the  run 

Integer  Set  5  - 

First  card  contains  501,  500  +  NTMS,  and  1  in  15  fields 
The  next  card(s)  contain  MISS(l)  through  MISS(NTMS)  in  15  fields 
12  entries  per  card  where  MISS(i)  is  the  mission  type  flown  on  the 
flights.>  NFLTS(i)  and  <  NFLTS  (i+l) 

Integer  Set  6  - 

First  card  contains  901,  900  +  NINSP,  and  1  in  15  fields 
The  next  card(s)  contain* s)  NIN(l)  through  NIN(NINSP)  in  15  fields 
12  entries  per  card  where  NIN(i)  is  the  number  of  flights  just 
preceeding  the  ith  inspection, 

integer  Set  7  - 

First  card  contains  921,  920  +  NRBWK,  and  1  in  15  fields 
The  next  card(s)  contain(s)  NRE(l)  through  NIN(NREWK)  in  15  fields 
12  entries  per  card  where  NRE(i)  is  the  number  of  flights  just 
preceeding  the  ith  rework 

The  next  cards  read  in  the  single  tables  as  required. 

TF  Table  - 

First  card  contains  1  and  NSTAB1  in  an  15  and  115  field  respectively 
where  NSTABt  is  the  number  of  points  in  the  TF  table, 

The  next  oard(s)  contain! s)  the  flight  number  entries  in  E10.3  fields 
six  to  the  card 

The  next  card(s)  contain! s)  the  times  corresponding  to  the  flight  number 
entries  in  E10.3  fields  six  to  the  card 
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PCB  Table  - 

First  card  contains  2  and  NSTAB2  in  an  15  and  115  field  respectively 

where  NSTAB2  is  the  number  of  points  in  the  PCB  table.  Note  that  if 

fodc 

the  PCB  table  is  derived  by  taking  n  points  from  the  function  p£ 
then  the  PCB  table  contains  n  points 

The  next  card(s)  oontain(s)  the  crack  length  entries  in  E10.3  fields 

J  Xj_,  x2i  *  *  *  ,  xn  are  abscissas  of  the 

(arranged  such  that  xi  +  1  >xi)  then  the  crack 

length  entries  are  x. ,  x  ,  ,<(1  x  . 

»  “  **  bcic 

The  next  cards  contain  the  ordinates  from  the  pr  function  in  E10.3 


six  entries  to  the  card, 
points  taken  from  p^e 


fields  six  entries  to  the  card, 
bde,  .  bde,  .  bde 
Pp  ) »  Pp  1  *  *  M  Pf 


Those  entries  are  ordered  as  follows: 


PI  Table  - 

First  card  contains  3  and  NSTAB3  in  an  15  and  115  field  respectively 
where  NSTAB3  is  the  number  of  points  in  the  PI  table  and  also  is  the 
number  of  points  taken  from  the  p^  function. 

The  next  oard(s)  contain  the  abscissas  of  the  points  taken  from  P° 
arranged  such  that  x., +1 >  x^.  These  numbers  are  input  in  E10.3 
fields  six  entries  to  the  card. 

L. 

The  next  oard(s)  oontain(s)  the  ordinates  of  tho  points  corresponding 
to  the  ahscissas  entries  above.  The  ordinates  are  entored  in  E10.3 
fields  six  entries  per  card. 


3CB  Table  - 

First  card  contains  4  and  NSTAB4  in  un  15  and  115  field  respectively 
where  NSTAB4  is  the  number  of  points  in  tho  SCB  table  and  also  is 
the  number  of  points  taken  from  the  function. 

U 

The  next  card(s)  contain(s)  the  abscissas  of  the  points  taken  from 
arranged  such  that  x^  +  1>  Xj .  These  numbers  are  input  in  E10.3 
fields  six  entries  to  the  card. 

The  next  card(s)  contain(s)  the  ordinutes  of  the  Sc  points  corresponding 
to  the  abscissa  entries  above.  The  ordinates  are  entered  in  E10.3 
fields  six  entries  to  the  card. 
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PCBM  Table  - 

The  PCBM  table  is  derived  within  the  Subroutine  PCBM  from  the  PCB 
table.  If  there  are  n  points  in  the  PCB  table,  then  there  are  2n-2 
points  in  the  PCBM  table.  The  first  point  of  PCBM  is  the  same  as  the 
first  point  PCB.  If  l<i<n  then;  the  abscissa  of  the  2i-2  point  of 
PCBM  is  the  same  as  the  abscissa  of  ith  point  of  PCB.  the  abscissa 
of  the  2i-l  point  of  PCBM  is  the  abscissa  of  ith  point  of  PCB  x  1.000001 
and  the  ordinate  of  the  2i-2  and  the  2i-l  point  of  PCBM  is  the  ordinate 
of  the  ith  point  of  PCB.  The  last  point  of  PCBM  is  the  same  as  the 
last  point  of  PCB. 


PDSBC  Table  - 

First  card  contains  1,  NDTAB11,  NDTAB12  in  an  15,  115,  and  110  field 
respectively  where  NDTAB11  is  the  number  of  stress  entries  in  the  table 
and  NDTAB12  is  the  number  of  different  missions  entered  in  the  table 
which  is  taken  from  the  Pg°  function. 

The  next  oard(s)  oontain(s)  the  stress  entries  x,, ,  x,,ni 
arranged  such  that  x^+1  >  x^  These  numbers  are  input  in  E10.3 
fields  six  entries  to  the  oard. 


The  next  card(s)  contain(s)  the  mission  number  entries  ,  ,y„,  »  , 

1  6 

yNDTAB12  arran8ed  auoh  yi+l>  yi’  These  numbers  are  input  in 

E10.3  fields  six  entries  to  the  card. 

t)C 

The  next  oards  contain  the  ordinates  of  the  Pg  function  arranged  so 
they  correspond  to  the  points  (x^y^),  (x^y^,'1',  (XNDTAB1 1 ,yl ^ ’ 


(xlty2) 


’  <XNDTAB11 '  1Wl)lZ)  f™“  th0  Cntr1'8  “b0Ve’  The‘e 


NDTAB11  •  NDTAB12  numbers  are  input  in  E10.3  fields  six  entries 


to  the  card. 


AB  Table  - 

First  oard  contains  2,  NDTAB21 ,  NDTAB22  in  an  15,  115  and  110  field 
respectively  where  NDTAB21  is  the  number  of  time  entries  in  the  table 
and  NDTAB22  is  the  number  of  different  missions  entered  in  the  table. 
The  next  card(s)  contain  the  time  entries  Xj,x2,  xNDTAB2i 

arranged  such  that  x^+1  >  Xj.  These  numbers  are  input  in  E10.3 
fields  six  entries  to  the  card. 
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WyiH! 


The  next  card(s)  contain(s)  the  mission  number  entries  ,  y^, 


yNDTAB22  arran8ed  such  that  "ivi  >  y^*  These  numbers  are  input 
in  E10.3  fields  six  entries  to  the  card. 

The  next  oard(s)  contain(s)  the  craok  lengths  arranged  so  that  thev 
correspond  to  the  pair  (x1,y1),  (xNDTAB2i*  > »  <x1»y2)' 

<XNDTAB21 '  yNDTAB22  ) ‘  These  NDTAB21  ’  NDTAB22  numbers  are  input  in 
E10.3  fields  six  entries  to  the  card. 


lit  i  a  .Mr. (aj, 


.a- tir  i’,ii  n  .V!  in'  nit 


5  EQUIVALENCE  TABLES 


The  technique  that  has  been  used  in  coding  this  routine  is 
to  place  all  input  and  output  numbers  in  blank  common.  All  input 
and  output  floating  point  numbers  are  called  parameters  and  are 
contained  in  P  (dimensioned  9000).  All  input  and  output  fixed  point 
numbers  are  oalled  integers  and  are  oontained  in  N  (dimensioned  1500). 
To  make  the  program  more  easily  read,  EQUIVALENCE  statements  are  used 
to  give  tne  P  and  N  number  more  recognizable  names.  The  RISKY  program 
parameter  and  integer  tables  are  given  below.  Also  given  below  is  a 
brief  description  of  the  tables  used  for  interpolation. 

Integer  Equivalence  Table 


N 

Dimension 

Term 

N 

Dimension 

Term 

1 

(1) 

IDENT 

51 

(5) 

NWT1(1) 

2 

(1) 

NPF1 

55 

NWT1 ( 5 ) 

3 

(1) 

NPF2 

57 

(2) 

NWT2 ( 1 ) 

4 

(1) 

NPF3 

58 

NWT2(3) 

5 

(1) 

NPF4 

6 

(1) 

NTH 

61 

(5) 

NTB1(1) 

7 

(1) 

NTW1 

65 

NTB1 ( 5 ) 

8 

(1) 

MONTH 

69 

(2) 

NTB 12(1) 

9 

(1) 

DAY 

70 

NTB12(3) 

10 

(1) 

YEAR 

72 

(2) 

NTB22( 1 ) 

11 

(1) 

NTI2 

73 
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12 
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NTW2 

75 

(15) 

NPW(l) 

13 
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NA 

89 

NPW (15) 

14 

(1) 

NCONTP 

15 
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NP 

16 
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NTMS 

101 
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17 

(1) 

NINSP 
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920 

NIN ( 20) 

22 
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NPRNT 

921 

(20) 

NRE(l) 

940 

NRE ( 20 ) 

961 

(400) 

NUMF(l) 

1360 

N'JMF(400) 

43 

(1) 

TTNSP 

44 

(1) 

JREWK 

45 

(1) 

KT 

46 

47 

(1) 

NP5IL 

48 

(1) 

NZERC 

49 

50 

(1) 

(1) 

NPAGE 

NDcrt: 

Parameter  Equivalence  Table 


Dimension 

Term 

U) 

DELCRK 

U) 

DELST 

(1) 

TIMAX 

(1) 

AI 

(1) 

ARR 

(1) 

ARRO 

P 

Dimension 

Term 

101 

(50) 

PSIL(l) 

ISO 

PSIL(50) 

151 

(400) 

PFBCT(l) 

550 

PFBCT(400) 

551 

(400) 

TIMCLC ( 1 ) 

950 

TIMCLC (400 ) 

951 

(20) 

TIMI(l) 

970 

TIMI (20) 

971 

(20) 

TIMR(l) 

990 

TIMR(20) 

991 

(100) 

PCBT(l) 

1090 

PCBT(IOO) 

1091 

(100) 

CRACKL ( 1 ) 

1190 

CRACKL( 100) 

1191 

(400) 

TIME ( 1 ) 

1590 

TIME (400) 

1591 

(400,15) 

PFBT( 1,1) 

7590 

PFBT(400, 15) 

7591 

(400) 

PFT(l) 

7990 

PFT(400) 

7991 

(400) 

PF(  1 ) 

8390 

PF(400) 

8391 

(400) 

PFF(l) 

8790 

PFF(400) 

8791 

(100) 

DAREAP(l) 

8890 

DAREAP( 100) 

50 

(1) 

RMISS 

51 

(1) 

PFBC 

52 

(1) 

DLAREA 
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Interpolation  Table  Listing 


Table 

Dimension 

Name 

Description 

TABL1(1,1) 
TABLl (200,1) 

(200) 

TF 

Flight  time  vs  number  of 
flights 

TABL1 (1,2) 
TABLl(200, 2) 

(200) 

PCB 

Crack  length  probability 
density  function  vs 
craok  length 

TABU  (1,3) 
TABLl (200, 3) 

(200) 

PI 

Inspection  function  vs 
crack  length 

TABL1(1,4) 
TABLl (200, 4) 

(200) 

SCB 

Critioal  crack  length 

function  vs  oraok  length 

TABLl (1,6) 
TABLl (200,6) 

(200) 

PCBM 

Modified  PCB  Table 

TABL2 (1.1) 
TABL2( 1000,1) 

(1000) 

PDSBC 

Stress  probability  density 
funotion  vs  stress  and 
mission  type 

TABL2 (1,2) 
TABL2 ( 1000 , 2 ) 

(1000) 

AB 

Crack  size  vs  time  and 
mission  type 

6  SAMPLE  PROBLEM 


The  following  sample  problem  is  presented  to  acquaint  the  user 
with  the  input  data  and  the  output  to  be  expeoted.  The  problem  has  a 
single  inspection  at  2200  flights  and  a  structural  rework  at  2500  flights 
The  input  oards  are  as  follows: 

207  1  0  7  0  4  1  8  18  1979  2  1  120  1 

1  30  1  1  1  50  68  0 

RISK  ANALYSIS  FOR  SEVERE  TRAINING  INCLUDES  INSPECT  AND  REWORK 
SINGLE  0.871  HOUR  MISSION  CRACK  SIZE  IN  INCHES 
1  6  1 

0.008  500.0  2265.0  0.080  0.030  0.040 

51  55  1 

11111 
87  58  1 

1  1 

75  78  1 

320 

101  130  1 

0  100  200  300  400  500  600  700  800  900  1000  1100 
1200  1300  1400  1500  1600  1700  1800  1900  2000  2100  2200  2300 
2400  2500  2600  2700  2800  2900 
501  530  1 

111111111111 
111111111111 
111  111 
901  902  1 

2200  3000 

921  922  1 

2500  3000 

1  2 

0.0  100.0 

0.0  87.1 

2  35 

33 
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0.0000316 

0.0000398 

0.0000501 

0.0000631 

0.0000794 

0.0001000 

0.0001259 

0.0001585 

0.0001995 

0.0002512 

0.0003162 

0.0003981 

0.0005012 

0.0005310 

0.0007943 

0.001000 

0.0012589 

0.0015949 

0.0019952 

0.0025119 

0.0031623 

0.0039811 

0.0050119 

0.0063096 

0.0079433 

0.0100000 

0.0125893 

0.0158489 

0.0199526 

0.0251189 

0.0316228 

0.0398107 

0.0501107 

0.0630957 

0.0794328 

0.0102 

0.0649 

0.361 

1.622 

6.385 

21.405 

61.154 

148.888 

308.86 

546.00 

822.51  1055.85 

1154.99  1076.62 

855.19 

578.86 

333.89 

164.11 

68.737 

24.533 

7.462 

1.934 

0.4271 

0.08037 

0.01289 

0.001762 

0.0002051 

0.2035-04 

0.1722-05 

0.1241-06 

0.7618-08 

0.3986-09 

0.1777-10 

0.0 

0.0 

3 

14 

0.00 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

0.10 

0.11 

0.12 

0.14 

1.00 

0.96 

0.85 

0.52 

0.36 

0.25 

0.19 

0.15 

0.11 

0.08 

0.06 

0.05 

0.04 

0.04 

4 

23 

0.001 

0.004 

0.005 

0.010 

0.015 

0.020 

0.030 

0.040 

0.050 

0.080 

0.100 

0.110 

0.120 

0.140 

0.160 

0.180 

0.200 

0.202 

0.204 

0.206 

0. 206 

0.210 

1.00 

200000.0  103000.0  95000.0  92000.0  75000.0  70000.0 

62000.0  57000.0  56000.0  61000.0  63000.0  70000.0 

66000.0  53000.0  43000.0  33000.0  25000.0  24000.0 

22900.0  21300.0  20000.0  0.0 


0.0 


1 


28 


1 


15000.0 

15700.0 

20000.0 

20500.0 

21000.0 

21500.0 

22000.0 

22500.0 

23000.0 

23500.0 

24000.0 

24500.0 

25000.0 

25500.0 

26000.0 

26500.0 

27000.0 

27500.0 

28000.0 

28500.0 

29000.0 

29500.0 

30000.0 

30500.0 

31000.0 

31600.0 

32000.0 

32500.0 

1.0 

0.0 

0.0 

807 . 2-06 

487.8-06 

365.9-06 

243.9-06 

130.66-06 

95.82-06 

34.84-06 

40.08-06 

33.10-06 

13.94-06 

12.18-06 

8.72-06 

15.68-06 

11.84-06 

8.36-06 

6.80-06 

2.96-06 

2.28-06 

1.38-06 

0. 878-06 

0.664  06 

i  0.348-06 

0.262-06  0.140-06 

i  0.112-06 

i  0.088-06 

2 

16 

1 

0.0 

10000.0 

12000.0 

14000.0 

14500.0 

15000.0 

15500.0 

16000.0 

16500.0 

17000.0 

17125.0 

17250.0 

17375.0 

17400.0 

17450.0 

18000.0 

1.0 

0.0001 

0.0010 

0.0015 

0.0025 

0.0030 

0.0060 

C .0140 

0.0260 

0.0520 

0.0860 

0.0960 

0.1110 

0.1280 

0.1500 

0.200 

0.205 

Based  on  this  input  the  following  ovitput  data  was  obtained. 
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APPENDIX  -  RISKY  PROGRAM  LISTING 


The  listing  given  below  is  a  FORTRAN  extended  language  routine. 
Seotion  3.3  gives  the  flow  diagram  for  the  various  subroutines. 


SUBROUTINE  INPUT  IKPT 

SUBROUTINE  FOR  INPUT  OF  FLOATING  MO  FIXED  POINT  NUMBERS  IKPT 

COMMON  P(9GBB)v  N (1508) ,  TTBLH29G.5I,  TAELZ(lBSa>2l  IKPT 
INTEGER  DAT,  TEAR  INPT 

DIMENSION  MTB116),  NT3I2t3J,  NTB22C31  IHPT 
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